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Matsumoto and Yoshimura have recently argued that the number density 
of heavy particles in a thermal bath is not necessarily Boltzmann-suppressed 
for T <C M, as power law corrections may emerge at higher orders in perturba- 
tion theory. This fact might have important implications on the determination 



of WIMP relic densities. On the other hand, the definition of number densities 



in a interacting theory is not a straightforward procedure. It usually requires 
renormalization of composite operators and operator mixing, which obscure 
the physical interpretation of the computed thermal average. We propose a 
new definition for the thermal average of a composite operator, which does 
not require any new renormalization counterterm and is thus free from such 
ambiguities. Applying this definition to the annihilation model of Matsumoto 
and Yoshimura we find that it gives number densities which are Boltzmann- 
suppressed at any order in perturbation theory. We discuss also heavy particles 
which are unstable already at T = 0, showing that power law corrections do 
in general emerge in this case. 
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1 Introduction 



The computation of number density of cosmological relics, such as weakly interacting 
heavy particles (WIMPs), neutrinos etc is usually realized by using a thermally average 
Boltzmann equation ffl, |2J. The scenario which emerges is that of a sudden freeze-out: 
the particle number density follows the equilibrium value until the freeze-out temperature 
below which the annihilation is frozen. The freeze-out temperature of pair annihilation 
model may roughly be estimated by equating the annihilation rate to the Hubble rate, 
the relic abundance is then computed by the thermal number density at that freeze-out 
temperature, thus, we expect it is Boltzmann suppressed if the freeze-out temperature 
(Tf) is much smaller than the mass of the particle. This is the typical case for a WIMP: 
the freeze-out temperature is typically 4-5 % of the mass of the particle namely, they are 
non relativistic at decoupling, and their annihilation cross section is of the right order to 
give a contribution to Qm of order unity, so that they are very good candidates for cold 
dark matter models ||. 

But from recent papers [Q, Matsumoto and Yoshimura (MY) have challenged the above 
conventional conclusion, claiming that two-loop corrections to the number density exhibit 
only power-law suppression in Tf/M, thus dominating over the Boltzmann-suppressed 
tree-level contribution if Tf is low enough. If confirmed, such a finding would imply that 
the present constraints on dark matter models are actually underestimated and should 
be carefully reconsidered. To be definite, MY introduced a toy model of two real scalar 
fields with Lagrangian 

±4 
4! 



-|x 4 -^V, (i) 

with M ^> m. The following relations were assumed among the coupling constants 

(A^^A 2 , |A|<|A X |<1, 

so that the light % particles act as an efficient heat bath for the heavy 0's. They then 
considered the quadratic part of the Hamiltonian for 0, 

H 4> = !0 2 + I(V0) 2 + iM 2 2 , (2) 

and, for temperatures T such that m ^ T ^ M they defined the number density of <fi 
particles as 

(£[<(, + counterterms)T 

N* = m (3) 

where MY define the thermal average for an operator A as 

= trAe^r _ (Om 
{ )T ~ tre- H / T (0|0> ' { ' 
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H being the total Hamiltonian. At tree-level, the well known Boltzmann-suppressed 
contribution is obtained 



N, 



(o) 



( 



MT\ 3/2 



V 2vr 



') 



-M/T 



The power-law contribution arises at two-loops and is given by 



N 



(2) 



y6 



69120 



(5) 



While in the standard case the contribution to Qm depends on the mass of the particle 
only trough logarithmic terms, in this case, with the different power law behavior fl5|), a 
simple calculation shows that the energy density goes like ~ M 1//4 , implying much stronger 
constraints on the mass of the particle. The interpretation of (|3|) as the physical number 
density of heavy (f> particles has been questioned by several authors |5|, ||, 0]: what emerges 
from all these papers is that the real problem is the definition of the number density for 
a heavy particle in a interacting theory. We will show that it is possible to give a non- 
ambiguous definition of the number density of an interacting particle at any temperature 
and at any order in perturbation theory. In the case of the MY model of eq. ([I]), our 
definition exhibits the expected Boltzmann-suppression at low temperature. We will find 
it very convenient to work in the real-time formalism of finite temperature field theory 
(RTF) ||, which has the advantage of keeping T = and finite temperature effects well 
separated. Indeed, as we will discuss, the ambiguities in the conventional definitions of 
the number density are mostly due to the T = O-renormalization procedure. Being able 
to single out T = ultraviolet (UV) divergences turns then out to be a great advantage 
when one is interested in genuine thermal effects. The tree-level propagator for a real 
scalar particle in the RTF has the structure: 



D = D^ + = zP[A ] + i(A - AS)A/"(M)B , (6) 

eynman propaj 

(exp(x/T) — (+)1)~ 1 and the two 



where A * } = (q 2 — M 2 + ( — )ie) 1 is the tree-level Feynman propagator, M the Bose- 
Einstein (Fermi-Dirac) distribution function, M{x) 
matrices P and B are defined as 



P[a(q)} 



B 



a(q) 

(a - a*)(q)9(q ) 

1 1 
1 1 



(a-a*)(q)e(-q Q ) 
-a*(q) 



(7) 



with 6{x) the Heaviside's step function. 
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The full propagator has the same structure as the tree- level one, eq. (0), with Ag 
replaced by its all-order counterpart 

II(go, \ q\) being the (renormalized) full self-energy at finite temperature. In the following, 
we will refer to zP[Ao] and zP[A] as the "T = 0" parts of the tree- level and full propa- 
gators, and to z(A - Ao)7V(|g |)B and z(A - A*)7V(M)B as the "thermal" parts. It is 
important to notice that what we call the "T = 0" part of the full propagator actually 
contains thermal corrections in the self-energy, so it is not the full T = propagator. 
The report presented here is a summary of our recent work ||. 



2 Renormalization of 

Definition of number density generally comes out from bilinear operators (see @ for MY 
model ), thus, we have to study composite operator renormalization. They arise when we 
consider local product of fields making up the operator. New UV divergences are induced 
when it is inserted in a Green function, which are generally not canceled by the Lagrangian 
counterterms. The renormalization of these divergences then requires the introduction of 
new counterterms. Renormalized operators can be defined, which are generically express- 
ible as linear combinations of all the bare operators of equal or lower canonical dimension- 
ality (for a thorough discussion of composite operator renormalization, see for instance 



TI|). Before considering the complete energy-momentum tensor, the analysis of the 
composite operator 0q (where the 0-label indicates the bare fields, masses and coupling 
constants) will be useful in illustrating the main points of renormalization and thermal 
averaging. We start by considering the T = case. In the model (|l|) the renormalization 
of <pQ then induces a mixing with Xo, an d the two renormalized operators, [<fi 2 ] and [x 2 ] 
can be expressed as 



[0 2 ] \ _ / l-5Z M 
[X 2 ] ~\ 




(9) 



In minimal-subtraction scheme (MS) the bare parameters can be expressed in terms of 
the renormalized ones as 

Ko-H A ^ 1 + ^ (D _ 4)fc ) > 

», 2 *rif, 6M2 \ *rif, ^ b k({*jh D )\ 




a -m'|l + ^-)=m-|l + ^g ^ } 'f ) 'l, (10) 



m 2 I \ m 2 
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where {A;} = {A^, A x , A}, D is the space-time dimensionality, and the important point is 
that the a^'s, 6^'s and 6™'s are mass, momentum, and temperature independent. Using 
eqs. (|I0| ) it is then possible to prove that, for any renormalized green function G the 
following relation holds 



M, 



dG 



2 dG 



dG 
dM 2 



+ mf 



dG 
dm 2 



showing that, being the right hand side manifestly finite, the insertion of the combination 
Mq0q + ttIqXo does not induce any UV divergence. The 5Z t 
constants can then be expressed in terms of the mass counterterms as 



SZ XX renormalization 



sz, 

5Z 



SM 2 
5m 2 



xx 



— ~ dZ <f>x 



(11) 



whereas the remaining counterterms 5Z^ X , SZ X ^, can be computed as 



5Z, 



4>x 



SZ- 



x<t> 



(xo( x )x(y)x(z)) 

{Xl{x)<j>{y)<t>{z)) 



DIV 



DIV 



3 Thermal average without mixing 

The mixing between the <fi 2 and x 2 operator makes the definition of a pure 0-number 
density quite cumbersome. A diagrammatic analysis can help us in identifying the T = 
nature of the pollution by the x field so that we can give a definition of thermal expecta- 
tion value for </> 2 which is free from it. Let us consider tr e~ H ^ T [(j) 2 ]/ tre~~ H / T . At tree-level 
in the RTF it is given by the two diagrams in fig. |l|a where the continuous line represents 
the T = part of the <fi propagator, the barred continuous line the thermal part and the 
double line the sum of the two. The first diagram is quadratically divergent and requires 
a subtraction, whereas the second one, thanks to the Bose-Einstein function contained in 
T>o(q), is finite. At 0(A) the dia grams of fig. |l|b contribute. There are also contributions 
O(A^) (such as higher order contributions containing also A x ) but their consideration is 
not necessary for our discussion, so we will limit ourselves to contributions containing only 
powers of A. The (quadratic) divergence due to the T = %-tadpole diagram (dashed 
line) is canceled by the mass counterterm. The only remaining divergence comes from the 
first diagram, in which both the lines, entering the cross are of T = type. It is the 
cancellation of this divergence which calls into play the mixing, via the renormalization 
constant SZ^ appearing in the last diagram. 
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(c) 



Figure 1: Diagrammatic expansion of ([4> 2 ])t- 



At 0(A 2 ) (fig. |l|c) the same happens, where: 



is the self energy to 0(A 2 ) renormalized by Lagrangian counterterms. The last diagram, 
containing a tree-level \ 2 operator multiplied by a 0(A 2 ) renormalization constant is re- 
quired in order to cancel the divergence coming from the integration over the momentum 
of particles flowing to the cross with T = propagators (the first diagram in fig. lc). 
Note that, due to the general relation 



i,i=l,2 ^ij 



0, Ejj 



being the self-energy with 
external thermal indices i and j, the contribution with both lines of thermal type also 
vanishes, once the summation over thermal indices is taken. The remaining loop diver- 
gences are canceled by usual Lagrangian counterterms and therefore require no mixing. 
The only divergence which is not canceled by neither Lagrangian nor composite operators 
counterterms is the tree-level one. Then, one usually defines the thermal expectation value 
of <j) 2 as the expectation value of [(f) 2 ] subtracted according to one of the two following 
ways @, 0, 



([4> 2 })t = 



tre~ H / T 



(O|[0 2 ]|O) 
(0|0) 



or 



tr [02] e -H/T ( O |[0 2 ]|O) 



tre-^/ T (0|0) 



(12) 

{ree fields 



The definitions above are both finite, and thus a plausible choice. However, both of them 
include the contributions from the first and the last graphs in figs. 0b and |l|c, that is, 
they require mixing. 

The above analysis suggests a new definition of thermal average of 2 , which represents 
the main point of this paper, namely 



<r<1 i(A - A*)(q)Af(\q \) , (13) 



(2?r) 4 

where D T (q) is any of the four component of the thermal part of the full propagator. 

Expanding A*-*) as A^ = A *^ + Aq IlAo + • • •, and recalling the relations between 
n(g , |gl) and S( ?0 , §, 

ReII( g ) = Re hnU(q) = e{q ) Im \E u (q) + E 12 {q)} , 

one can check that the perturbative expansion of fig. [I] is reproduced, apart from the 
'problematic' diagrams. The definition in ([13|) then automatically gets rid of all the 
unpleasant features of the conventional definitions in eq. ([TJ), namely, the need of an 
arbitrary subtraction and, more importantly, composite operator renormalization and 
mixing. 



4 Averaging the energy- momentum tensor 



11, 12 . Ex- 



The renormalization of the energy-momentum tensor is discussed in refs. 
pressing it in terms of bare parameters, the operator 

V = PfaVfa + d"xod u Xo ~ + • • ■ , 

is finite when inserted in a Green function, that is, it does not require extra counterterms 
besides those already present in the Lagrangian. The dots in the formula above represent 
pole terms proportional to (d^d" — g^d 2 )^ and to (d^d u — g^d^xl- Since we are 
interested in the thermal average of T^ v alone, translation invariance insures that such 
pole terms do not contribute. On the other hand, a subtraction of the T = contribution 
as in eqs. fll2j) is needed. Working with renormalized parameters, T^ v is expressed as 

TV = T C Z + T& , (14) 

where Tj£L has the canonical form and T^ v t contains only the Lagrangian counterterms. 
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The finiteness of T^ v means that the divergences induced by, e.g., the composite op- 
erator 4> 2 x 2 are canceled by the Lagrangian counterterms and the particular combination 
of operators 4 , x 4 > {d<f)) 2 , ■ ■ ., appearing in T^ v ' . In general, if we split T^ u in parts, each 
part separately will require composite operator renormalization. Indeed, this is what is 
done by MY in ref. H. They split the total Hamiltonian as 



T 



oo 



tot 



H<t> + H x + H int + H c , t , , 
and compute thermal average of (given in eq. (2)). In the RTF it gives 
tr^e"^ T , f d 4 q 



tYe -H/T 



2q 2 Q - (q 2 - M 2 ) 



(2vr) 4 
[D T=0 + D T 



12 



(?)] 



(15) 



where, differently from eq. (|13j), also the T = part of the full propagator appears. It 
is important to recall here what we have noticed after eq.(H), namely that the T = 
full propagator contains thermal effects via the self-energy. Such contributions are not 
canceled by purely T = subtractions like those defined in eq. dl2|) . By perturbatively 
expanding [D T=0 ]i2 in ( |15"D at 0(A 2 ) we find, apart from divergent contributions which 
require composite operator counterterms to be added to H^, the "famous" power-law 
contribution 

d 4 q 



i(A, 
= div 



(2^)4 
d 4 q 

(2vr) 4 

2 




M 



T=Oi 



12 



2g 2 - (a 2 

2q 2 Q - (q 2 - M 2 
A* 2 )Ren-(A 2 + A* 2 )e(g )Imn 



o(A 2 ) 



-Qo) 



(16) 



1 A 2 T 6 



+ 



69120 M 2 

where e(q ) = 9(q ) — 9(—q ). The same contribution, with opposite sign, is found from 
the corresponding piece of the thermal average for H x . 

Paralleling our discussion in the previous section, we will instead define the energy 
density of the 4> field as 



d 4 q 



M 



[D- 



12 



(?) 



(17) 



and analogously for p' . The complete energy density will then be split as 



Ptot = {HtotjT 



tiH tot e- H l T (0\H tot \0) 
tre~ H / T (0|0) 

= P's + Py + P' 



int ) 
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k 

Figure 2: Two-loop contribution to the self energy . 

so that the contribution of eq. fll6|) enters p[ nt instead of p'^. Notice that neither p'± nor 
p' x are affected by the T = subtraction. 

The definition in ( |18D exhibits three remarkable properties: 

i) The splitting in ( |I8J ) is closed under operator mixing and composite operator renor- 
malization. No composite operator counterterm is required to make p 1 ^ or p' x separately 
finite and, since p to t is also finite, the same is true also for p[ nt . Thus, our definition of 
number density, differently from that considered by MY, is independent on the renormal- 
ization scale p []. There is no "pollution" from p' x to p'^ at any value of p. 

ii) None of the three pieces in (|T8|) contains 0(A 2 T 6 /M 2 ) terms. Indeed, the contri- 
bution of eq. (0) and the analogous one for the go both into p' int and there, because 
of their opposite sign, cancel. 

in) p'j, is Bolt zmann- suppressed for T <C M. Indeed, writing [D T ]i 2 (g) explicitly 

[D ]M = - ( g2 -M 2 -Ren( 9 )) 2 + Imn(g)^ (l90|) ' (19) 

we see that, due to the Bose-Einstein function, non-Boltzmann-suppressed contributions 
to the integral in eq. (|17]) might only come from momenta |go| ^ T (^C M). But, for 
these values of momenta, it is the imaginary part of the self-energy (Figure 0) which is 
Boltzmann-suppressed. It can be easily understood looking at Figure 3 : 

q ^ ' • v q 

° - ° (a) 




Figure 3: a) Decay, b) Annihilation, c) Landau damping . 

for |go| ^ T(^C M) only the annihilation and the Landau damping contribute. In the 

*Of course, at any finite order in perturbation theory we still have the usual /i dependence induced 
by Lagrangian counterterms. 
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former case (fig. |3|.b) the on-shell particle in the initial state comes from the heat-bath 
and then carries a M factor. In the Landau damping case (fig. [|c), the energy to create 
the on-shell <p has to be provided by the \ from the heat bath, so that a Boltzmann- 
suppressed Af has to be payed in this case too. 



5 A Better Definition for N„ 



Using the 0-energy density as given in eq. fllTD , we can now define the number density 
for T C M in a way analogous to MY's definition, eq. (§), that is 

N <t> = ^ for T<M. 

However a better definition, valid at any value of T, and again based on the average of a 
quadratic operator, can be given as 

N * = I ^M^Uiq) ■ ( 2 °) 

As p^, eq. (p0|) exhibits the nice features of not requiring composite operator renormal- 
ization and of being Boltzmann-suppressed at low temperatures, but in addition it may 
be extended to high temperatures or, equivalently, to the massless limit. The origin of 
eq. (|20|) can be derived in analogy to the case of a complex scalar field, $, in case an exact 
U(l) symmetry is imposed to the theory. The number density for the real scalar field 
should differ from the previous quantity in two respects. First, particle and antiparticle 
contributions should be summed up, in order to account for their 'Majorana' nature. This 



is obtained by taking the modulus of \q Q \ inside (|20|). Second, a factor 1/2 is needed, in 



order to take into account the number of degrees of freedom. 



6 The Decaying Case. 

The discussion of the previous sections is straightforwardly generalizable to any model in 
which heavy particles annihilate into lighter ones, both bosons and fermions. The basic 
conclusions remain unalterated, the main reason being the Boltzmann-suppression of the 
imaginary part of the self-energy at small momenta. 

The situation changes drastically if the heavy particle is unstable at T = 0. If, for 
definiteness, it decays into two lighter particles of mass m < M/2 the propagator exhibits 
an imaginary part for 2m < po < M already at T = 0, so that there is no reason to expect 
that it is Boltzmann-suppressed for T < M. The heavy particle is actually a resonance, 
whose energy-momentum can take values much lower than the peak at p 2 = M 2 , of course 
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paying the usual Breit-Wigner suppression far away from the peak. In this case, power- 
law contributions to the resonance number density do in general emerge, as we will briefly 
discuss. 



We will consider for definiteness the annihilation model discussed in ref. [13] and in 
the first of refs. [|J , in which the interaction between the heavy $ and the light x~ bosons 
is due to the term 

C int = -^ X \ (21) 
where we assume that fi <C M. In this model, the imaginary part of the self-energy is 



oc \x , so that, using (|l^) in (|20|) we get 



N, ~ ^T 5 , (22) 



with T = fi 2 /(32irM) the T = decay rate on-shell. In ref. [13| a different behavior 



~ T 4 was claimed. Indeed, the following argument clarifies the physical origin of the 
power-law contribution and confirms the T 5 behavior. 

In the thermal bath, unstable $ particles are continuously produced in x~X annihila- 
tions. Being the $ a resonance, the production energy, u, can be much smaller than M. 
The produced $'s eventually decay with inverse lifetime T = fi 2 /(32nuj) = YqM/uj. Thus, 
the number density of $'s of energy around u, n*(u;), obeys the rate equation 

— j — Tn (u), (23) 
where 7 is the rate per unit volume of the process 

XX ~~ * & ~ > everything . (24) 
The above equation leads to the equilibrium value 

$/ \ 7 

n eq{v) = f ■ 

The production rate is given by 7 ~ N^av, where the inclusive cross section a may be 

A 

computed using the optical theorem as a ~ /i 2 ImA$/s, s being the square of the $ 
four-momentum, and A$ the full $ propagator. 

Taking m \fs ~ T M, the initial x particles are relativistic (i.e. N x ~ T 3 ) and 
we get a ~ n 4 /(M 4 T 2 ). Putting all together, we get a contribution to the total $-number 
density from the region u ~ T of order 



<^- T ^w T ^w^ b - (25) 



The $'s with energies closer to the peak of the resonance live longer but their contribution 
to the number density is Boltzmann-suppressed, as two highly energetic x' s are required 
in the initial state to produce them. 
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As the energy of the relevant $'s is typically 0(T), the energy density is 
p$ ~ O {j^3^ 6 ^j + 'Boltzmann — suppressed! , 

and p$/p x vanishes as T 2 at low temperatures. Before concluding this section, we observe 
that the definition of a sensible number density along the lines of eq. (§), i.e. starting 
from the free Hamiltonian, does not seem quite appropriate in this case. Indeed, as we 
have just discussed, the width of the $ plays a crucial role, allowing energy values much 
lower than the tree-level mass-shell. It seems that, if one wants to insist in looking for 
definitions like eq. (|]), a free Hamiltonian for quasi-particles, along the lines discussed for 
instance in ref. jnj, should rather be used. 



7 Conclusions 

The definition of particle number density in a interacting theory is a delicate matter. The 
necessity of giving a meaning to divergent composite operator calls into play operator 
mixing, so that a clear separation between different particle species turns out to be a 
renormalization scale dependent procedure. We showed that in the RTF, it is possible to 
give a proper definition of number density free from this problem. In the framework of 
MY annihilation model , our definition of the 0-energy density is Boltzmann suppressed, 
no power suppressed terms of the type found by MY appear. This implies that the compu- 
tations of number density of cosmological relics, on which dark matter models are based, 
are correct. 

For an unstable particle , we showed, that power law contribution to the equilibrium num- 
ber density do in general emerge and this may have important cosmological consequences 
for example in relation to the generation of the Baryon Asymmetry. It is clear that in 
this case the usual Boltzmann equation have to be modify to include off-shell effects. On 
this subject and the possible cosmological consequences of this new scenario, we are still 
working on [O . 
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